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ABSTRACT 

Using mean field MHD, we discuss the energetics of optically thin, two tempera- 
ture, advection-dominated accretion flows (ADAFs). If the magnetic field is tangled 
and roughly isotropic, flux freezing is insufficient to maintain the field in equipar- 
tition with the gas. In this case, we expect a fraction of the energy generated by 
shear in the flow to be used to build up the magnetic field strength as the gas flows 
in; the remaining energy heats the particles. We argue that strictly equipartition 
magnetic fields are incompatible with a priori reasonable levels of particle heating; 
instead, the plasma @ in ADAFs (= gas pressure divided by magnetic/turbulent 
pressure) is likely to be ^ 5; correspondingly, the viscosity parameter a is likely to 
be <, 0.2. 

Subject headings: accretion, accretion disks - hydromagnetics - turbulence 



1. Introduction 

Recent work (Narayan & Yi 1994, 1995ab; Abramowicz et al. 1995; Chen et al. 1995; 
Nakamura et al. 1997; Manmoto, Mineshige, & Kusunose 1997; see Narayan, Mahadevan, & 
Quataert 1998b, or Kato, Fukue, & Mineshige 1998 for reviews) has revived interest in a class of 
hot optically thin accretion solutions first discovered by Ichimaru (1977): advection-dominated 
accretion flows (ADAFs). In ADAFs, the accreting gas is unable to cool efficiently and most 
of the energy generated by turbulent stresses is advected onto the central object. As a result, 
the gas heats up to nearly virial temperatures and adopts a two-temperature configuration 
(Shapiro, Lightman, & Eardley 1976; Rees et al. 1982), with the ions significantly hotter than 
the radiating electrons. 

In this paper, we derive and discuss the energy equations for particles, magnetic fields, and 
turbulence in a two-temperature ADAF (§2). There are several mutually exclusive prescriptions 
in the literature for treating the energetics of ADAFs (compare, e.g., Esin, McClintock, & 
Narayan 1997, hereafter EMN and Manmoto et al. 1997). The differences between these 
prescriptions lead to significant differences in the predicted spectra of accreting black holes 
(compare, e.g., Narayan et al. 1998a and Manmoto et al. 1997 for Sgr A*), so it is important 
to resolve these discrepancies. This is undertaken in Appendix A. 
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Non-relativistic ADAFs exhibit a curious tension between two conflicting requirements. 
On the one hand, energy advection causes the flow to be virial, with a profile of density and 
temperature (T oc p 2//3 ) that makes a non-relativistic gas nearly isentropic. On the other hand, 
dissipation of turbulent energy should cause the gas entropy to increase inwards. This tension 
is particularly relevant for optically thin, two temperature, ADAFs which have seen the most 
success in applications to real systems (see Narayan et al. 1998b for a review); for such flows, 
Ti 3> T e and the gas is quite accurately approximated as non-relativistic. 

In an effort to alleviate this difficulty, Narayan & Yi (1995b; NY95) suggested that a 
proper treatment of the flow dynamics should explicitly (if phenomenologically) account for 
the turbulence/magnetic fields in the flow (since these components are not expected to behave 
like non-relativistic particles). Using the energy equations derived in §2, we explore this sug- 
gestion in some detail. In §3 we show that the most straightforward implementation of NY95's 
suggestion (a constant "/?" for the flow, which is utilized by many workers in the field) is still 
incompatible with the turbulent heating of particles. In §4 we suggest an alternative method of 
including turbulence in the dynamics of ADAFs which explicitly accounts for particle heating; 
we present analytical and numerical solutions consistent with this scenario. Finally, in §5 we 
summarize and discuss the implications of this work. 



2. The Energetics of ADAFs 



The goal of this section is to derive energy equations for particles, magnetic fields, and 
turbulence in a two-temperature ADAF (see §18 of Kato et al. (1998) for a related derivation). 
Later in this paper (and in Appendix A) we use these equations to clarify several issues relevant 
to spectral and dynamical models of ADAFs. 

Strictly speaking, since ADAFs are nearly collisionless, we should use the Boltzmann 
equation to analyze their dynamics /energetics. We assume, however, that the use of MHD is 
justified by the presence of the (turbulent) magnetic field. In particular, this assumes that 
various wave-particle interactions (pitch angle scattering) maintain approximate isotropy in 
the distribution functions. With this simplification, the equations governing the structure of 
the accretion flow are 

^ + V-( P v) = 0, (1) 

/°isr + p( v ■ v ) v = ~ VP + ffi + r( VxB ) xB + ^V 2 v, (2) 
at 47r 

9B 

— = V x (v x B) + 7/ B V 2 B, (3) 
at 

prf = H - q -, (4) 

where p, v, P, g, B, s, rj v , and tjb represent the mass density, velocity, gas pressure, gravitational 
acceleration, magnetic field strength, particle entropy per unit mass, viscous diffusivity, and 
magnetic diffusivity in the flow (respectively). In equations dU)-(@), d/dt denotes a Lagrangian 
derivative while d/dt denotes an Eulerian derivative. Equation (||) is the sum of the electron 
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and ion energy equations, H is a volumetric particle heating rate (due to viscous, resistive, or 
plasma dissipation mechanisms) and q~ is the volumetric cooling rate. Below we will discuss 
the individual ion and electron energy equations. 

We emphasize that H is the heating rate of the particles. In ADAFs, H is not necessarily 
equal to q + , the rate at which energy is generated by turbulent stresses/shear in the accretion 
flow, since some of q + may be used to build up the magnetic field strength and turbulent kinetic 
energy as the gas flows in (see §2.3). Similarly, s in equation (|j) is not the total entropy of the 
fluid, but only that of the particles. 

We use mean field theory to derive energy equations for the turbulent flow (see, e.g., 
Parker (1979) or Speziale (1991) for general discussions of mean field theory and turbulent 
closure schemes; Balbus & Hawley (1998) and Kato et al. (1998) give nice applications of this 
formalism to accretion theory). We assume that each quantity characterizing the flow can be 
expressed as a sum of an ensemble- averaged ("mean") quantity and a fluctuating quantity, e.g., 
v = U + u, B = Bo + b, P = Pq + p, where the fluctuating quantities (u, b, p) have zero mean 
(we denote ensemble-averages by ()). We assume that the turbulence is roughly incompressible, 
but we do not take V-v = 0. Rather, we make the minimal modification to the incompressibility 
assumption which allows us to treat the compressive heating of particles/field/turbulence; we 
thus take a static density field, p = p(r) (i.e., we take dtp = 0, but not dtp = 0). Another 
simplification of this work is that we take Bo = 0, i.e., we assume that there is no mean 
magnetic field. 

ADAFs are a one scale problem; the vertical scale height (H v ) is of order the local radius 
in the flow and the inflow time of the gas is of order the rotational period (S7™ 1 ). Assuming 
that turbulent eddies have sizes ~ H v and turnover times ~ Q^ 1 , this makes the use of a mean 
field theory somewhat suspect. For now, however, it is the best we can do. 

By subtracting the mean versions of equations (|2|) and (||) from the full equations, one 
derives equations for Ui and 6j, the fluctuating velocity and magnetic field. Multiplying these 
equations by Uj and bj, respectively, yields equations for the evolution of the Maxwell and 
Reynold's stress tensors. The resulting energy equations for the fluctuating magnetic field and 
turbulent kinetic energy are given by (see also §18 of Kato et al. 1998) 

d t (b 2 ) = -2(b 2 )d k U k + 2(b i b k )d k U l 
- d k [(u k b 2 ) - T] B d k {b 2 )^ 

+ 2{bib k d k Ui) - (b 2 8 k u k ) -2r jB ((d k b i ) 2 ) (5) 

and 

d t (pu 2 } = -{pu 2 )d k U k -2{puiU k )d k Ui 

- 2d k {^-{pu 2 u k ) + (p tot u k ) - (^jty - |^ 2 }) 

- 2 (j^ d kUi) + 2{p tot d k u k ) - 2r] v ((d k Vi) 2 ), (6) 

where p to t = b 2 /8-K +p is the total fluctuating pressure and d t = d t + U k d k is the Lagrangian 
derivative with respect to the mean flow. The physical meaning of these terms is as follows. 
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The first line on the right hand side of equations (||) and (^) represents the generation of 
energy via compression and coupling to shear in the mean flow. The second line gives the 
divergence of the flux of energy due to turbulence and microscopic resistivity/viscosity. The 
three terms on the third line represent exchange of energy between the velocity and magnetic 
components of the turbulence, fluctuating PdV work due to correlations (if present) between 
the fluctuating pressure and the divergence of the turbulent velocity field, and dissipation by 
microscopic resistivity /viscosity, respectively. 

Taking the mean of equation (||) yields the mean energy equation for the particles, 

dt(e) = -j(e)d k U k + (H)-(q-) 

- d k (^^j-{pd k u k ), (7) 

where e = P/{"f — 1) is the thermal energy per unit volume (and hence (e) = -Po/(7 — 1)) and 
7 is the adiabatic index of the electron + ion gas; if both species are non-relativistic, 7 = 5/3, 
if both are relativistic, 7 = 4/3, etc.Q 

A number of different "7's" will be used in this paper: 7 refers to the adiabatic index 
of the electron + ion gas, while j e and 7, refer to the individual electron and ion adiabatic 
indices, respectively. The quantity jt refers to the "adiabatic index" of the turbulence, which 
is a measure of how efficient flux freezing/compression is at increasing the turbulent energy 
(§4). Finally, 7 9 refers to the total adiabatic index of the fluid, containing contributions from 
the particles and turbulence (§2.2). 



2.1. The Total Energy Equation 



An equation describing the total energetics of the magnetic field, turbulence, and particles 
is obtained by summing equations ©-(f?]), which yields 

d t E + d k F k = G - (q-) (8) 

where 

£ = (^ 2 ) + (£) + <£> (9) 

is the total (turbulent plus thermal) energy per unit volume, 

/ u k b 2 \ 1 9 7 
F k = {-^) + 2(puu k ) + —^(puk) 

' —d k {u z ) - —dk(b z ) (10) 



4vr / 2 8vr 



1 In equation (fjj) we have written the particle entropy gradient as (pTdts) = d t (e) — j(e)Ukdk In p. We use 
this interchangeably with (pTd t s) = pd t (cf/(7 — 1)) — (ct)UkdkP, where c? a = P/p is the sound speed. 
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is a turbulent flux (the viscous/resistive contribution is negligible) 



G 




TikdkUi 



(11) 



is the term associated with generating energy, and 



T ik = (puiU k ) 



47T 



(12) 



is the stress tensor. 



2.2. Approximations 



Equations (|5|)-(|T2]) are fairly rigorous energy equations for the particles/fluctuations. They 



are, however, too complicated for some of the calculations of interest to us (e.g., calculating 
the dynamics of an accretion flow). We therefore make a number of approximations to the 
mean field energy equations which make them more suitable for simple analytical or numerical 
calculation. These simplified equations maintain (hopefully) a number of the essential features 
of the mean field equations. 

First, we must provide a closure relation for the third order turbulent quantities (i.e., the 
turbulent flux). We choose, without much justification, to ignore these terms, setting = 0. 
Second, we assume that the diagonal components of the turbulent stress tensor are equal, in the 
sense that {b^} ps (bg) ~ with analogous equations for {vf). We do not, however, ignore 
the off diagonal terms in the stress tensor since they are responsible for angular momentum 
transport. The generalization to arbitrary diagonal components for the stress tensor is given 
in Appendix B. 

We also assume that the flow is azimuthally symmetric, has Uq = 0, and rotates uniformly 
on spherical shells (i.e., Uj, = rf2(r) sin where the rotation rate Q is independent of 9; see 
Narayan & Yi 1995a). Taking Uq = requires that the angular scale height of the flow be 
roughly independent of radius (Narayan & Yi 1995a; Abramowicz et al. 1997). With these 
approximations, we can write, in spherical coordinates, 



This equation demonstrates that, under compression, the velocity component of the turbulence 
behaves like a gas of adiabatic index 5/3 while the magnetic field behaves like a gas of adiabatic 
index 4/3. This is a consequence of our assumption of isotropic turbulence. 

We express the level of turbulence in the flow in terms of the gas pressure via two param- 
eters: 




(13) 



0b 




and /?, 




(14) 



(b 2 /8ir) 
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We generally take (3b = (3 V = (3 (i.e., the turbulence is "Alfvenic"). Our definition of (3 is that 
utilized in the plasma physics literature. A number of workers in the accretion literature (e.g., 
NY95) define a u f3" via /3 a d v = Pgas/Ptoti i- e -> the fraction of the total pressure supplied by 
the gas. This is related to our (3 by /3 a d v = 3/3/(3/3 + 1) or /3 a d v = /3/(/3 + 1), depending on 
whether one defines the magnetic pressure to be b 2 /2A-K or b 2 /8n. 

Finally, we set 6 = vr/2, i.e., midplane values, in the above expressions; this is a crude - 
but standard - form of height integration. With these assumptions it is rather straightforward 
to show that the total energy equation for the flow (eq. ||]) takes the form (dropping the () 
for simplicity) 

pU r d r ( — — J = U r vfd r p + q + - q~ , (15) 

where vt, the effective isothermal sound speed, is related to the real isothermal sound speed of 
the gas (c s ) and the turbulent sound speed (q) by 

, f 2 = ^ + c? = ^(l + /r 1 ), c* = ^, c? = |^ = |. (16) 

The quantity 7 5 is an effective adiabatic index for the flow which includes contributions from 
the particles, the magnetic field, and the turbulence. For isotropic Alfvenic turbulence, we 
have 

_ 7 + 3 (7-i)r 1 , 17 , 

^~ i + 2 ( 7 -l)/3-l- [U) 



If we take the ion + electron gas to be nonrelativistic (a good approximation) , equation (|17|) 
simplifies to 

5/3 + 6 

* - wn- < 18 > 

This expression for j g is slightly different from that of Esin (1997) who found, in our notation, 
7 5 = (15/3 + 8)/(9/3 + 6). This is because she neglected the ram pressure of the turbulence, 
i.e., she set (3 V — > oo. 

The particular expressions for vt and 7 g in equations (16) and ( |l7| ) are a consequence 
of our assumption of isotropic Alfvenic turbulence, the validity of which is highly uncertain. 
Appendix B generalizes these expressions to the case of arbitrary diagonal components for the 
stress tensor. 

The quantity q + in equation (jl5|) is the rate of energy generation in the flow due to 
coupling between the turbulence and the shear. It is given by 

q+ = ~ Tr * r lP (19) 
which can be "closed" using a Shakura & Sunyaev (1973)-like relation, 

T r<t> = -apv t H v r—. (20) 

Equation ( |l5|) demonstrates an important point. In order to write an energy equation 
with q + (the standard viscous heating term) as a source term, one must consider all energy 
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components in the flow (hence the appearance of the "effective" sound speed and "effective" 
adiabatic index in eq. [jig] , which contain contributions from the turbulence and the particles). 
It is, strictly speaking, incorrect to write a particle energy equation with q + as a source term, 
since the particle heating rate, H, is, in general, somewhat different from q + . The relationship 
between the "particle" and "total" energy equations is discussed further in Appendix A. 



2.3. The Electron, Ion, and Turbulence Energy Equations 

The mean electron and ion energy equations are given by the individual versions of equa- 
tion (0) which, utilizing the assumptions given above (steady state, azimuthal symmetry, drop- 
ping third order turbulent quantities, etc.), simplify to 

U r d r ei = 'jieiUrdr In p + Hi - q ie - g~ (21) 

and 

U r d r e e = j e e e U r d r In/? + H e + qie - 9^, (22) 

where H e (Hi) is the energy dissipated out of the turbulence which heats the electrons (ions), 
7 e (7i) is the monatomic ideal gas adiabatic index for the electrons (ions), q~ (q~) is the 
electron (ion) cooling rate, and qi e is the ion-electron energy exchange rate. We have dropped 
the () for conciseness. 

As it stands, however, equations ( |2~I|) and (|22| ) are not closed since there is no prescription 
for calculating H e and Hi. Let us define 5h to be the fraction of the dissipated turbulent 
energy which heats the electrons, so that H e = 5hH and Hi = (1 — 5h)H. To get a handle on 
the particle heating rate, H, we note that it is equal to the total energy dissipated out of the 
turbulence (magnetic and velocity), and thus appears in the energy equation for the turbulence. 
Summing equations (||) and (|6|), and employing the various simplifying assumptions outlined 
in §2.2, we find that the energy equation for the turbulence simplifies to 

UrSr g + ^ = (|£ + Urdr ta , + ,+ _ H , m 

which becomes 

2pU r d r c 2 t = c 2 t U r d rP + q + — H (24) 

for isotropic Alfvenic turbulence (for which b 2 /8ir = pu 2 /2, c 2 = b 2 /8irp, and the turbulent 
"adiabatic index" is 1.5). 

Equations (|23| ) and ( |24"1) state that the rate of change of the turbulent energy as the gas 
accretes inward is that supplied by compression/flux freezing (the d r p term) and coupling to 
the shear in the mean flow (g + ) minus the energy dissipated into particle heat (H). We can 
rewrite these relations schematically as 



Qt = q + - H, 



(25) 
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where Qt is the rate of change of the "turbulent entropy." Note that H, the particle heating 
rate, shows up as a cooling term in the turbulent energy equation. 

In ADAFs, just as one cannot ignore the advection term in the particle energy equation, 
one cannot set H = q + and Q t = in the turbulent equation. We do, however, expect that 
H ~ q + (rather than, say, H <C q + or H = q + ). This is because of the one scale nature of 
the turbulent energy equation: the timescale on which shear generates energy (the rotational 
period) is also the timescale on which energy cascades to small scales to heat particles (i.e., 
the eddy turnover time), which is also ~ the timescale on which equipartition attempts to be 
established (the inflow time of the gas).0 

For the particular case of isotropic turbulence, Q t ^ and thus H q + . This is because, 
under compression, isotropic Alfvenic turbulence behaves like a 7 = 1.5 gas, i.e., b 2 oc p 1 ' 5 
(cf eq. p4| ). By contrast, for a quasi-spherical flow, equipartition with the gas requires 
b 2 ~ pc 2 oc p 5 / 3 (§3). Compression (flux freezing) therefore does not supply sufficient energy 
to keep the turbulence in equipartition with the gas. Some of the energy generated by the 
coupling between the turbulence and the mean flow (q + ) goes into increasing the turbulent 
energy, while the remainder heats the particles (hence H 5a q + ). 

For non- isotropic turbulence configurations, one could have H <^ q + . This occurs when 
the turbulence behaves like a gas of adiabatic index J> 5/3 (see Appendix B). Flux freez- 
ing/compression then acts to increase the level of turbulence above the equipartition value. In 
this case, the turbulent energy dissipated in the flow (H) must exceed that generated by shear 
(q + ), in order for the flow to remain gas pressure dominated (Bisnovatyi-Kogan & Lovelace 
1997). 

Many models of two-temperature ADAFs in the literature parameterize the turbulence 
by assuming a constant (3 for the flow (note that in deriving the above equations we have not 
employed this assumption). In this case, one can substitute c 2 = c 2 j3~ l into equation (|24j) and, 
knowing the dynamics of the flow (which gives c 2 , p, U r , q + , etc., as a function of r), calculate 
the particle heating rate that is implicitly assumed by these models (since H is the remaining 



unknown in equations [23] and ]24|). This is carried out in the next section; surprisingly, and 
problematically, it turns out that many models in the literature require H <C q + , i.e., most of 
the energy generated by turbulent stresses (q + ) is implicitly assumed to be advected radially 
by the turbulence (Qt), rather than heating the particles (H). 



3. The Dynamics of ADAFs 



In this section we reanalyze current treatments for phenomenologically including turbu- 
lence/magnetic fields in the dynamics of ADAFs. We show that explicitly decomposing the 



2 By contrast, in thin accretion disks the inflow time of the gas is much longer than the the rotational period. 
To an excellent approximation, then, one can set Qt = and H — q + . This is the turbulent analog of being 
able to drop the entropy gradient in the particle energy equation, setting q~ = H (which becomes the more 
standard q~ = q + by the above argument). 
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energetics of the flow into a turbulent/magnetic component and a gas component, as done in 
the previous section, leads to a significant reinterpretation of these models. 

Consider a steady axisymmetric accretion flow. The dynamics of viscous hydrodynamic 
accretion can (following Narayan, Kato, & Honma 1997, hereafter NKH) be described by the 
following height integrated equations representing conservation of mass, radial momentum, 
angular momentum, and energy, 



±(pr 2 H U r ) = O, (26) 
ar 

U r ^-^r = -nlr--U P vf), (27) 
dr p dr 

Ur d(nr 2 ) _ 1 d f apvyH e dn\ 
r dr pr 2 Hg dr \ Qk dr j 

di 

pU r T-=q + -q-. (29) 
dr 

Anticipating the quasi-sphericity of the flows of interest, we have carried out height integration 
on spherical shells, so that Hg is the angular scale height of the flow; where necessary, we take 
H = vt/(rSl K )fl 

In the spirit of the previous section, we interpret the pressure (sound speed) and entropy in 
equations (^6|)-(p9|) as effective quantities, containing contributions from both the particles and 
the turbulence. In particular, we take v 2 = c 2 s + cf as the sound speed (eq. [pp and interpret 
equation (29) as equation (|l5|) (with the effective adiabatic index j g ). It is straightforward 



to employ the averaging/approximation scheme discussed in the previous section to the MHD 
momentum equation (eq. Q) and derive equations (p7|) and (|28|) . 



3.1. Self-Similarity 



In the limit of q~ C g + , cooling provides a rather minor perturbation when calculating 
the dynamics of the flow. One can therefore set q~ = (1 — f)q + in the energy equation (with 
/ ~ 1) so that fq + of the viscous energy is advected by the flow. In this limit, NY94 showed 
that, if 7 5 and / are independent of radius, and Qk is Newtonian, equations (f26"|)-(f2S|) admit a 
self-similar solution, with the various flow variables being power laws in radius (note that this 
entails that (3 is a constant, independent of radius). In particular, Hq = constant, U r oc ar -1 / 2 , 



3 In spherical coordinates, balancing pressure and the centrifugal force in the direction implies p oc 
exp [— 0.5fi 2 r 2 cos 2 #/w 2 ] (e.g., Narayan, Barret, & McClintock 1997), assuming v t and to be independent of 9 
(this, but not Q,, vt independent of z, is a reasonable approximation for ADAFs). Defining He to be the spherical 
average of the density, i.e., Hq = (47r) _1 § sm 6 d9d<j>p yields Hg = (tt/2) 1 ^ 2 x erf (x -1 /\/2), where x — vt/(yir). 
This is the correct expression for He for ADAFs. It is, however, numerically cumbersome, so for now we stick 
with the (more standard) expression given in the text. 
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v 2 oc r 1 , p oc a r' " 3 / 2 , and Q oc (5/3 - -/g) 1 / 2 Q K (see NY94 or Narayan et al. 1998b for the 
complete scalings). 

Of particular interest for our purposes is the dependence of on j g . Non-relativistic 
(ig = 5/3) advection-dominated accretion flows tend to be non-rotating. The reason is that 
the above equations describing the structure of the accretion flow in ID exhibit the following 
"singularity." Approximate hydrostatic equilibrium in the radial direction implies a virial 
flow, for which v\ ~ r 2 ^ 2 ^ oc r _1 . Together with angular momentum conservation, this implies 
U r ~ or -1 / 2 . Mass conservation then implies p oc a _1 r~ 3 / 2 . Together, then, mass, radial 
momentum, and angular momentum conservation impose a T oc v 2 oc p 2 ^ density/temperature 
profile on the flow. The remaining flow variable, f2, adjusts itself to ensure that the energy 
equation is satisfied. The profile T oc p 2 ^ is, however, an adiabat for a non-relativistic flow 
with 7 9 = 5/3. This means that there can be no entropy generation as the gas flows inwards, 
which requires that f2 and q + go to zero. 

The remainder of this paper is predicated on the assumption that this singularity as 
7 5 — > 5/3 is, at some level, real, and is not an artifact of the approximate treatment of 
the dynamics given above. Height integration schemes differing from ours do not show this 
behavior (e.g., Chen et al. 1997; Manmoto et al. 1997). We believe that these models, based on 
cylindrical averaging, are not accurate representations of the quasi-spherical flows of interest, 
but only 2D or 3D calculations will fully resolve this issue. 

NY95 proposed that non-relativistic ADAFs (such as the much applied T{ T e variety) 
behave differently from a 5/3 gas because of turbulence (in particular, magnetic fields) in the 
flow. The "formal" justification for this suggestion is given by the derivation of §2 (cf eqs. fl5| - 
[ |lg| ] and the accompanying discussion). In what follows, we employ the most straightforward 
implementation of NY95's suggestion; we analyze the dynamics of ADAFs assuming non- 
relativistic particles and a constant (3 turbulent /magnetic contribution to the pressure and 
energy density of the flow. 

In this case, the effective adiabatic index 7 9 ^ 5/3, and the flow can be differentially rotat- 
ing. This is important since we expect real flows to possess non-negligible angular momentum. 
If, however, we take the self-similar solution for j g ^ 5/3 and compute the entropy gradient for 
a nonrelativistic gas (7 = 5/3), we find Tds/dr\ gas = 0. Accounting for turbulence/magnetic 
fields to make the flow as a whole (7^) behave differently from a 5/3 gas does not alter the 
fact that the gas component of the flow still has 7 ~ 5/3 (see Appendix A), and therefore has 
no entropy gradient in the self-similar regime. It therefore does not alleviate the fundamental 
tension between being virial and the need for particle heating. 

Using the energy equations given in §2 (see the discussion below eq. |24||), it is straightfor- 
ward to show that H = q~ = (1 — f)q + and Qt = fq + - This implies that, although including 
constant turbulence in the dynamics of a non-relativistic ADAF allows for differential rota- 
tion and entropy generation (q + 7^ 0), it does so by requiring that almost none of the energy 
generated by turbulent stresses heats the particles (since for / ~ 1, H <C q + ). Instead, this 
energy is assumed to be stored as turbulent energy in the flow, i.e., it is advected by the tur- 
bulence, rather than the particles. For the self-similar solution, this conclusion is independent 
of /?. For larger (3, the turbulent energy is smaller; 7 9 is, however, closer to 5/3 and so Q and 
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q + are smaller as well. The solution always adjusts itself so that the energy generated by the 
shear is precisely what is required to maintain constant (3 in the flow. 

The condition H <^ q + means that only a small fraction of the turbulent energy cascades 
to small scales and heats the particles. This seems unreasonable since intuition and dimensional 
analysis dictate that ADAFs should have H ~ q + (although not H = q + ; see §2.3). 



4. An Alternative Dynamics: Explicit Inclusion of Particle Heating 



One way to resolve the above contradiction is to give up the assumption of constant /3. 
Instead, let us assume that the magnitude of the turbulence/magnetic fields in the flow adjusts 
so as to allow the requisite level of particle heating. 

To analyze this scenario, we replace the single energy equation representing "particles plus 
turbulence" (eq. |]l5|) with the following two energy equations, the first for the particles and 
the second for the turbulence, 

pU r d r (^rrl = U r (? s d r p + H-q~ (30) 

pU r d r [^jj = U r c 2 t d r p + q + -H. (31) 
The quantity 7 (74) is the particle (turbulent) adiabatic index. For c 2 = b 2 /8irp and jt = 1-5, 



equation (31) was derived in §2.3 (it is equivalent to eq. ^4|); here we allow for a slightly more 
general fluid model of the turbulence by considering a general value of 74 • Equation ( |30|) is the 
sum of the individual particle energy equations (eqs. |2l| and p2f). 

To ensure particle heating, we replace the assumption of a constant (3 with the assumption 
that a constant fraction 77 of the energy generated by turbulent stresses heats the particles, 
i.e., 

7] = H/q + = constant. (32) 

(3{r) is now an output of the dynamical equations. Mass, radial momentum, and angular 
momentum conservation are still given by equations (|26|)-(p8|), taking the total pressure to 
be the sum of the gas and turbulent pressures, v 2 = c 2 + c 2 . To be consistent with previous 
treatments, we define q~~ = (1 — f)q + (although q~ = (1 — f)H might be a more physical 
paramatrization) . With these definitions, 

H - q- = (77 + / - l)q + and q + - H = (1 - rj)q + . (33) 

As discussed in Appendix B, a crucial assumption of the present analysis is that j t ^ 5/3. 
This is equivalent to the statement that the effect of compression/flux freezing on the turbulent 
energy is insufficient to maintain the turbulence in equipartition with the gas. In this case, it 
is q + , the coupling between the turbulence and the shear, which is primarily responsible for 
increasing the turbulent/magnetic energy as the gas accretes inwards. 
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4.1. Self-Similarity 



We search for analytical solutions to equations (|2q)-(|28|) and (|30|)-(|3l[), assuming that /, 
r], and a are independent of r. For simplicity we take a 2 <C 1 (almost certainly a reasonable 
approximation), f3 ^ 1 (also likely to be a reasonable approximation), and 7 = 5/3 (non- 
relativistic gas, i.e., ion dominated). The assumption on (3 is needed in order to express all 
quantities in the solution as power laws in radius (in particular, it allows us to approximate 
c 2 = i^/(l + /? -1 ) as c 2 « v 2 (1 — (3~ 1 ))- With these assumptions it is relatively straightforward 
to obtain the following approximate solution to the above equations. At a fiducial radius ro, 
let P = 0q. For r < ro, we then have 

vt = ^rQ K , U r = ~(3/2-a)rn Ki (34) 

/ r \ o-3/2 / \ -2a 

n = n °{-J , = #>(-) • ( ; «) 




.5/% (3/2 

2(7 7 + /-l)/3 + (l-7 7 )( 7t -l), (37) 
(5-3 7t )(l-/ 



6( 7f -l)(77-l)-4/' 



(38) 



and 



„ = (l-/)(5-3 7f -4a)+6 a ( 7 «-l) , . 

I 5_3 74 _ a(4 _ 6(7t _i )) ' ^ ) 

For a = 0, 77 = 1 — /, and the above solution reduces to that of NY94. Equations (|34|)-(f3S|) are 
therefore the generalization of the NY94 self-similar solution to the case when particle heating 
is explicitly included in a flow with non-relativistic particles. We note that, even though the 
flow has vl oc /5 2//3 , which was the origin of the difficulty in NY94's self-similar solution, particle 
heating is accounted for in the present solution; this is because (3 7^ constant, so that neither 
the gas nor the turbulence is strictly virial (although their sum is). 

The key feature of the above solution is that (3 increases as a function of decreasing r (a > 
0). In the self-similar solution of NY94 all of the energy generated by turbulent stresses goes 
into maintaining the turbulence in equipartition with the gas (regardless of the magnitude of 
(3). Explicit inclusion of particle heating therefore must (and does) yield subthermal magnetic 
field strengths. Once some of the energy goes into particle heating, there is insufficient energy 
available to maintain constant (3 (given the constraints imposed by being virial). The increase 
in (3 is accompanied by a decrease in O/f^ (recall that the NY94 solution has constant 
f2/f2#). The total pressure and radial velocity (and hence density), however, have the same 
(virial) radial dependence as in NY94. 

A priori, we expect H ~ q + . We therefore take r\ = 1/2 (equal turbulent and particle 
heating) as a canonical value. For / = 1, this yields a = 1/14 and (3 oc r -1 / 7 for 7i = 1.5 and 
a = 1/6 and (3 oc r _1//3 for jt = 4/3 (isotropic magnetically dominated turbulence, for which 
c t = b 2 /24irp). For an ADAF with a radial extent ~ 10 4 Schwarzschild radii, (3 changes (over 
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the whole flow) by a factor ~ 4 (~ 20) for jt = 1-5 (4/3). Explicit inclusion of particle heating 
can therefore lead to noticeably subthermal magnetic fields. The precise numerical factor by 
which (3 changes is, however, quite sensitive to the (poorly understood) choice of jt an d rj. 

There are two unsatisfying features of this analytical solution. The first is the strong 
dependence on the outer boundary condition; the radial power law for (3 implies that we 
must specify its value at some fiducial radius; this in turn determines the value of (3 in the 
interior (where the emissivity peaks). Secondly, the analytical solution assumes that, for / = 1, 
constant (3 corresponds to rj = (no heating of particles). While correct for the self-similar 
solution of NY94, this is not strictly correct for global models of ADAFs (see below). This 
motivates us to consider global solutions of equations (p^)-(|28|) and (|3Ql)-(j3l"l). The unsettling 
features of the analytical solution are not reproduced by these numerical calculations, although 
the general inference of subthermal magnetic fields is. 



4.2. Global Dynamics 



We have calculated global (transonic) ADAF models by solving equations (pq)-(p8D and 
(|30|)-(|3l]) numerically using a relaxation method. Here we discuss the numerical method used; 
in §4.2.1 we give sample solutions. The implications of these calculations are discussed in §5. 
The models discussed in this section are the generalization of NKH's global calculations to the 
case when particle heating is explicitly included. 

We work in the pseudo-Newtonian potential introduced by Paczyhski & Wiita (1980), for 
which £l 2 K = GM/[r(r — r g ) 2 ], where r g is the Schwarzschild radius and M is the mass of the 
black hole. Integrating equation (^) implies 1 = —pU r r 2 HQ, where we normalize the density 
by taking the accretion rate to be M = 4-7T. Using equation (p6|) and the expression for Hq 
from above, we integrate equation (p8|) once to find 



dn u r n K (nr 2 -j)_ 



dR ar 2 v 2 

where j is the specific angular momentum accreted by the black hole 



Equations (p7|), (pC|)-(|3l|), and ( |40| ) are the four first order differential equations which we 
solve (mass conservation simply normalizes the density). We also have two eigenvalues, j and 
r s , the location of the sonic point. We therefore require 6 boundary conditions. 

Three of the boundary conditions are the values of c s , and q (or equivalently, Q, vt and 
(3) at the outer boundary (which we usually take to be at r out = 10 5 r 5 ). We typically specify 
(3 at the outer boundary and use the self-similar solution to calculate $7, c s , and q. 

Equations ©-(H) and ©-© can, using standard techniques, be manipulated to show 
the presence of a sonic point (at r = r s ), at which din \U r \/dr = N/D, where 



2p + 2 U 2 
p + 2 



D = ^i~^, (41) 
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y _rra-tf)_„A_.W + tf (42) 



dr / p£/ r t> 2 (p + 2)' 

/ = (7 - l)(r? + / - 1) + (7* - - 1), (43) 

and 

y _(7-lW + (-B-lW (44) 

A well-behaved solution must have N = D = at r = r s , which provides two additional 
boundary conditions. The location of the sonic point is not predetermined but is our second 
eigenvalue (it is typically at ~ 5r g for a ~ 0.3). 

The final boundary condition is the torque free condition, j = $7r 2 , at the horizon. For 
numerical reasons, we typically apply this condition at r = 1.5r g , but we have confirmed that 
the results are unchanged if we apply it at r = l.lr g . Applying this final boundary condition 
at the horizon is a defect of our acausal treatment of viscosity (see NKH for a discussion of 
this point and Gammie & Popham 1998 for a causal treatment of viscosity in ADAF models). 



4.2.1. Results 



For the numerical results in this section we set / = 1, 7 = 5/3, 74 = 1.5, and a = 0.3, and 
explore the solutions as a function of rj and the value of j3 at the outer boundary. We focus on 
understanding the evolution of (5 as a function of r, and the relationship between (5 and rj. 

To start with, we apply the self-similar boundary condition corresponding to f3 = 1 at 
f = r out . Figure la shows the radial velocity (U r ) and total sound speed (vt) for models 
ranging from 77 = to rj = 1. Figure lb shows the density for the same models. As predicted 
by the self-similar solution, the basic dynamical quantities (total pressure, density, etc.) are 
independent of the details of the particle heating (i.e., rj). The insensitivity of p, vt, and U r 
to rj is a generic feature of our models; we therefore do not show additional plots of these 
quantities. 

Figures lc and Id show Q/Qk and f3 as functions of r for rj varying from to 1 in steps of 
0.2. Here there is a nontrivial dependence on rj. In order to understand the behavior of f3(r), 
the solid lines in Figure 2a show f3(r) for an alternative outer boundary condition, /3 = 15 at 
r = r ou t (all other parameters are the same as in Figure 1). 

The self-similar solution of §4.1 predicts that j3 should be a monotonically increasing 
function of r, namely, /3 oc Pout r ~ 2a where, for a fixed 74, a > is just a function of rj (see 



eq. |38|). For rj ^ 1, this is not always seen in the global calculations. In fact, in contrast to 
the self-similar solution, (5 shows a tendency to approach a roughly constant value at small r 
(especially near the sonic point, inside of which the fluid flows in nearly adiabatically) , and to 
decrease with radius if it starts off with a large value (Fig. 2). 



Deviations from self-similarity in the density and temperature profiles of the flow, driven 
primarily by the existence of a sonic point, which modifies the radial velocity and therefore 



-15- 



p, can lead to more noticeable changes in the entropy gradient (because the entropy gradient 
involves a cancellation between the temperature and density gradients). Since the value of ij is 
oc the entropy gradient (cf eqs. Q|-[^2[), it, and hence 0, is sensitive to such deviations. This 
is the basic reason why self-similarity is less accurate when calculating r] and /?, even though 
it is reasonably good for p, cj, etc.. 

Nonetheless, the key qualitative conclusion of the self-similar solution is confirmed by the 
global calculations. Explicit inclusion of particle heating at reasonable levels (r) ~ 1/2) leads 
to subthermal magnetic fields. In fact, there is an approximate mapping between the value of 
(3 in the interior of the flow and 77, which is independent of the outer boundary condition (this 
is in contrast to the self-similar solution, which was a power law in r, and thus very sensitive 
to the outer boundary condition). This is seen in Figure 2b, which shows /?(r) for r\ = 1/2, 
taking (3 ou t = 1,5,15, and 50. In the inner portions of the flow, where the emissivity peaks, 
the solutions roughly converge to a common (3 ~ 5. 

Furthermore, the global calculations confirm that strong levels of turbulence in the interior 
of the flow are incompatible with reasonable levels of particle heating. Strictly equipartition 
fields {(3 = 1) near r ~ 1 are only obtained for unphysically small values of r/ ^ 0.1 (Fig. Id 
and 2a). 



5. Discussion 

In §2 of this paper we derived the energy equations for particles, magnetic fields, and tur- 
bulence in optically thin, two temperature ADAFs (using mean field MHD). The fundamental 
energy equations are given in equations @ - (|l2]); for explicit calculations, however, we have 
used drastically simplified versions which approximate the turbulence as an additional fluid 
in the problem (the approximations required are detailed in §2.2). One application of this 
analysis, given in Appendix A, is to resolve a disagreement in the literature between several 
treatments of the energetics of ADAFs. 

The primary focus of this paper, however, has been to clarify and explore a previously 
underappreciated piece of physics regarding optically thin, two-temperature ADAFs. The 
issues addressed stem from the fact that non-relativistic ADAFs are, at some level, a singular 
problem. Strong advection leads to a virial flow in which non-relativistic particles tend to be 
isentropic. This is because, at least in the self-similar regime, being virial requires T oc p 2 / 3 , 
which is an adiabat for non-relativistic particles with 7 = 5/3. On the other hand, the turbulent 
heating of particles should cause the entropy of the gas to increase inwards. To resolve these 
conflicting requirements, either the flow must be non-rotating (unlikely) or else a formalism 
must be developed which allows for both differential rotation and particle heating. 

This issue is particularly relevant for optically thin, two temperature, ADAFs which have 
been extensively applied to observed systems; for such flows, Tj » T e and the gas is quite accu- 
rately approximated as non-relativistic. Following NY95, a number of workers have suggested 
that, because of the "singularity" a proper treatment of the flow dynamics should include the 
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turbulent / magnetic contribution to the pressure and energy density (the electron contribution 
is expected to be less important if (3 ^ T{/T e , which is well satisfied in most models). The 
derivation in §2 "formalizes" this suggestion. 

Previous attempts at phenomenologically including turbulence in the dynamics of ADAFs 
assumed a constant (3 for the flow (e.g., NY95, EMN, Narayan et al. 1998a). What was not 
appreciated, however, is that this still does not account for the heating of particles; rather, 
the implicit assumption is that most of the energy generated by shear in the flow is stored 
as turbulent energy (§3). This energy is in fact precisely what is required to maintain the 
constant (3. 

The scenario required by these models is unlikely to be realized in real flows. Strongly 
turbulent plasmas inevitably heat particles (energy cascades to small scales, etc.). To assess 
the implications of this, we impose particle heating in the dynamical equations by specifying 
that some fraction (= r/) of the energy generated by turbulent stresses heats the particles. 
The remaining energy (a fraction 1 — rj) is used to build up the magnetic field strength and 
turbulent kinetic energy as the gas flows in. The parameter rj therefore replaces the parameter 
[3. The turbulence/magnetic field in the flow {(3) is now an output of the dynamical equations; 
it evolves so as to maintain the specified level of particle heating. 

For thin accretion disks, one expects 77 = 1 to good approximation due to the mismatch 
between the timescale on which the gas pressure increases (the inflow time) and the timescale 
on which shear acts to increase the turbulent energy (the rotational period). For a not <C 1, 
no such mismatch exists for ADAFs and so rj 7^ 1, although we do expect 77 ~ 1 (rather than 
v « i; §2-3). 

We have given analytical and numerical solutions for the dynamics of ADAFs at constant 77. 
These are, respectively, the generalization of the NY94 self-similar solution and NKH and Chen 
et al.'s (1997) global calculations (in the Paczyhski potential) to the case where particle heating 
is explicitly included. The results of the self-similar solution are particularly straightforward 
to understand. In previous, constant (3, models, the constancy of (3 was achieved at the 
expense of particle heating. If, instead, we require a reasonable level of particle heating, 
less energy is available for the turbulence and so (3 must increase as the gas flows in. The 
turbulence/magnetic field therefore become subthermal. We emphasize, however, that the 
predicted density, radial velocity, and total pressure of the flow are relatively unchanged (with 
respect to NY94 or NKH) by the present analysis (§4; Figure 1). 

Self-similarity requires that (3 monotonically increase as the gas flows in, while our global 
calculations show that this is not reproduced if the turbulence is somewhat subthermal. The 
global calculations do confirm, however, that strong levels of turbulence ((3 ~ 1) are incom- 
patible with reasonable levels of particle heating (since they require rj <C 1; see Figures Id and 
2a). If we prioritize particle heating, taking rj ~ 1/2, we find that (3 converges to ~ 5 in the 
interior of the flow, regardless of its initial value at large radii (Figure 2b). 

Explicitly calculating (3 requires understanding the non-linear saturation of the turbulence; 
this would provide the full profile, r/(r). Our point is that, regardless of the details of the non- 
linear saturation, one expects that 77 is not 1. Taking a lower bound of rj ~ 1/2 yields the 
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corresponding bound, (3 ^ 5. 

The (3 used in this paper refers to both the magnetic energy in the turbulence and the 



turbulent kinetic energy (cf eq. [14]); in particular, we assume roughly Alfvenic turbulence, for 
which these two energies are equal. Our analysis suggests that consistency with particle heating 
requires both subthermal magnetic and turbulent kinetic energies; on dimensional grounds, i.e., 
independent of the physics of angular momentum transport, this requires a relatively smaller a. 
This can be seen by noting that a « —T r( p/pc 2 s in a Shakura & Sunyaev-like closure relation, 
where c s is the sound speed. Prom the definition of the stress tensor (eq. |]l2fl ), we expect 
T r( j) pc 2 , where c 2 = c 2 //3 is the square of the turbulent velocity. Dimensional analysis 
therefore implies a 1//3. Quantitatively, then, the constraint ^ 5 suggests that a ^ 0.2. 



5.1. Caveats 



To conclude, we discuss several key assumptions implicit (or explicit) in our analysis. 

First, the correct value for rj(r) is highly uncertain, and the inferred (3 depends non- 
trivially on the rj chosen (see Figures 1 & 2). For example, one might think that 77 and (3 
would be anti-correlated (i.e., preferential build up of the turbulence if it is subthermal). In 
principle, one could solve the mean field energy equations (eqs. || - 0) using, e.g., a second 
order closure scheme, and thus calculate 77 (Speziale 1991), but this is well beyond the scope 
of our paper. 

Second, we treat the turbulence/magnetic fields in the flow using a crude fluid model. 
Even within the context of this model, we assume that the effect of flux freezing/compression 
on the turbulent energy is insufficient to bring the turbulence into equipartition with the gas 
(jt & 5/3). This requires that less than half of the magnetic energy is in radial perturbations 
and that more than a third of the turbulent kinetic energy is in radial perturbations (Appendix 
B) . Our physical picture is that the primary source of the turbulent energy is shear in the mean 
flow. 

A noticeable omission is our neglect of third order (and higher) turbulent quantities (§2.2). 
These represent, among other things, the turbulent diffusion of energy /entropy, and thus are 
quite relevant for the issues discussed in this paper. Unfortunately, turbulent transport is 
particularly complicated in ADAFs, since both convection and MHD turbulence are believed 
to be important. 

Although we consider a radially varying (3, we take a to be a constant. Aside from 
simplicity, our motivation for this restriction is that there is generally no self-similar advection- 
dominated accretion solution when a decreases as the gas flows in (which would be expected 
by our analysis if a oc [3~ l ). This is because a decreasing a leads to an increase in density 
more rapid than r~ 3 / 2 (if a oc r 9 , U r oc r _1 / 2+9 , and p oc r~ 3 / 2-9 ). The compressive heating 
associated with this density profile is sufficiently large that it generally yields a decreasing 
entropy gradient for the flow. The only viable self-similar solution in this case is one which is 
cooling dominated, rather than advection-dominated. 
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Finally, our calculations are (as is standard) a one dimensional (height integrated) analysis 
of a three dimensional problem. Relaxing this restriction will likely make the singularity 
inherent in the dynamical equations less severe. The extent to which this will alleviate our 
concerns about particle heating is, of course, unclear. 
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A. Related Work 



Here we compare our analysis of the energetics of ADAFs with several other treatments 
in the literature, and clarify several issues about the use of "effective" adiabatic indices. 

EMN treat the energetics of an ADAF using both a total energy equation and an electron 
energy equation. Our derivation/discussion of the total energy equation in §2 is a more rigorous 
treatment of the problem, which arrives at conclusions similar to theirs (which is based on Esin 
1997). The one difference is that they neglect the turbulent velocity (i.e., (3 V — ► oo), so that 
their expressions for c 2 and "f g are slightly different (the most general expression for 7 5 is given 
by eq. [B6] of Appendix B). 



EMN's treatment of the electron energy equation (their Appendix A) is, however, quite 
different from ours. They take the electron internal energy and pressure to be effective quan- 
tities, containing contributions from the magnetic field. They thus argue that the adiabatic 
index which shows up in the electron energy equation is also an effective 7 (in fact they take 
7 e = 7 fl ). The physical picture on which this is based is that (1) all of the energy generated by 
turbulent stresses heats the particles, i.e., H = q + , and (2) the energy required to bring the 
magnetic energy from its flux-frozen value to the equipartition value comes from the thermal 
energy of the particles (hence 7 e smaller than its monatomic value). The problem with this 
analysis is that it incorrectly treats the thermal energy of the particles as a source of magnetic 
energy. Within the context of MHD, turbulence and magnetic fields can act as a source of 
thermal energy, while the converse is not possible (this can be seen formally using the mean 
field energy equations in §2). If we set H e = (no turbulent heating of electrons), then the 
turbulence and magnetic field should not appear in the electron energy equation (except as a 
third order turbulent diffusion term), while in EMN's treatment they still would (through the 
effective 7). To conclude, then, the appropriate adiabatic index in the electron energy equation 
is not an "effective" one, but rather that of a monatomic ideal gas. 

Manmoto et al. (1997) and Nakamura et al. (1997) also solve the energy equations for 
ADAFs. They take electron and ion equations analogous to ( |2~l|) and (|22| ) as their basic energy 
equations, with H e = 0, Hi = q + , and 7 e = 7; = 5/3. Setting H e = assumes no turbulent 
heating of electrons, which we regard as overly optimistic. Setting 7 e = 5/3 is correct in spirit 
(in that there is no "turbulent" component to 7 e ), but it assumes non-relativistic electrons; 
this is incorrect in the inner ~ 10 2 — 10 3 Schwarzschild radii. Setting Hi = q + assumes that 
all of the energy generated by turbulent stresses heats the ions. As emphasized in this paper, 
however, H is, in general, 7^ q + (although it is likely ~ q + ). 

There is a close relationship between the "ion" and "total" energy equations discussed 
above and in §2. In particular, in treating the "ion/total" energy equation, one can either 
use a monoatomic adiabatic index or an effective adiabatic index. If the former, the heating 
term should be H 7^ q + , since one is truly writing a particle energy equation. If the latter, 
the heating term should be q + (see eq. ]l5| ) . Put another way, it is precisely the inequality of 
H and q + which allows one to use an effective adiabatic index. This can be seen explicitly by 
"solving" the turbulent energy equation (eq. |24|) for H and substituting it into the particle 
energy equation (eq. ) , in which case one readily derives the total energy equation (eq. [jl5| ) . 
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No such analysis is allowed for the electron energy equation, however, since it is a standard 
assumption of ADAF models that T e <C T% and the electrons couple poorly to the turbulence 
(Sh <C 1). The uncertainty introduced by H ^ q + is therefore absorbed into an uncertainty in 
(the already uncertain) electron heating parameter, 5h, rather than in the adiabatic index. 

The above arguments regarding the electron adiabatic index in ADAFs break down if 
the distribution functions are anisotropic. In this case adiabaticity cannot be described by a 
single adiabatic index (cf double adiabatic theory). In particular, for anisotropic distribution 
functions, thermal energy can be a source of magnetic energy (the fire-hose instability being a 
well known example), contrary to the arguments using MHD given above. ADAFs are nearly 
collisionless systems; isotropy is therefore not guaranteed. Nonetheless, we (as with all other 
workers in the field) have proceeded on the assumption that it is a reasonable approximation. 
Anisotropic distribution functions are known to be unstable to generating modes (e.g., Alfven 
waves, fast modes, whistlers) which act to restore isotropy (e.g., Melrose 1980). The timescale 
for such processes is typically of order the cyclotron period, much shorter than the dynamical 
time in the accretion flow. 



B. Anisotropic Turbulence 



When simplifying the mean field energy equations in §2, we assumed that the diagonal 
components of the turbulent stress tensor are all equal (i.e., ~ b 2 ^ ~ bg, etc.). This is equiv- 
alent to assuming that the "turbulent pressure" is isotropic. In this Appendix, we generalize 
the derivation of §2 to arbitrary diagonal components for the stress tensor. 

Assuming azimuthal symmetry, Ug = 0, and uniform rotation on spherical shells, but 



allowing for arbitrary values of T rr ,Tgg, and T^, the energy generation term of equation (11), 
which shows up in the total energy equation for the fluid, becomes 

G = - ((^ 2 ) + 2 (£) + ^ e >) dkUk ~ T ^ ~ ( Tee + ~ G ° d > ( B1 ) 

where 

dfl 

G d = sin Or— - Tg^Q cos 6 (B2) 

is the contribution to G from the off-diagonal components of the stress tensor. This is the 
same as derived in §2.2. 

For T rr = Tqq = Tgfa, as taken in §2, equation ( |Bl| ) reduces to equation (|lg|). This is 
because dU r /dr + 2U r /r = V-U under the assumptions of Ug = and azimuthal symmetry. In 
this case, the terms associated with geometrical convergence of the flow in spherical coordinates 
(the dU r /dr and U r /r terms) are equivalent to a term oc dp/dr, i.e., a PdV work term. 

As equation ( |B1| ) shows, however, one cannot in general write the "geometrical" terms as 
"compressive" terms and thus one cannot derive a fluid approximation for the turbulence in the 
flow. Assuming self-similarity, U r oc r~ 1//2 , allows one to do so. In this case, dU r /dr ~ — V-U/3 



and U r /r ph 2V • U/3. With this simplification, equation (Bl) reduces to 



- 21 - 



G 




dkUk - G d, 



where 



(V 2 r) 



and hh 



(b 2 r) 



(B3) 



(B4) 



(v 2 ) (b 2 ) 

are the fraction of the turbulent kinetic energy (h v ) and magnetic energy [hb) in radial per- 
turbations, respectively. 

Expressing the level of turbulence in the flow in terms of the /3's of equation (|14|) allows one 



to derive a total energy equation for the flow which is identical in form to equation (15). The 
expressions for the total sound speed, Vt, and the effective adiabatic index, j g , are, however, 
different. They are now given by 



and 



(B5) 



7 + ( 7 - 1) ftr 1 (I - 2h v ) + ( 7 - 1) /y 1 (I + 2h b ) 



(B6) 



l + (7-l)/3,- 1 + (7-l)/? b - 1 
If we set h v = hi, = 1/3 and j3\, = (3 V = j3, equations (|B5| ) and (|B6| ) reduce to equations (|l6| 
and (|l7|), respectively. 



Equation (B6) demonstrates the simple, but important, point that the effect of flux freez- 
ing/compression on the turbulent energy depends on the geometry of the turbulence. If we set 
in equation ( |B6| ) , we recover the "adiabatic index" associated with flux freezing of the 
magnetic field, namely 

(B7) 



lb 



2 

3 + 2h b . 



The corresponding adiabatic index for the velocity component of the turbulence is 

lv = 7 --2h v . (B8) 

It should be emphasized again that equations flB3|)-([B8|) are strictly valid only in the self-similar 
regime. 

In the bulk of this paper we have assumed that flux freezing is insufficient to keep the 
magnetic field/turbulence in equipartition with the gas, i.e., that 7^ ^ 5/3. By equation (|B7[) , 
this can be recast as an assumption on the magnetic field geometry, namely h b ^ 1/2, i.e., less 
than half of the magnetic energy is in radial perturbations. The corresponding constraint on 
the turbulent velocity field is that h v ^ 1/3, i.e., more than 1/3 of the turbulent kinetic energy 
is in radial perturbations. 

Both of these constraints (hb ^1/2 and h v <; 1/3) are satisfied in numerical simulations of 
thin accretion disks (cf Hawley, Gammie, & Balbus 1996). These are, however, of questionable 
relevance for the accretion flows at hand. For laminar, non-rotating, spherical accretion, it 
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is known that hb ~ 1: the magnetic field becomes purely radial (Shapiro 1973), and flux 
freezing is energetically important, invalidating the analysis of this paper. This result is also 
of questionable relevance for the accretion flows at hand since, while quasi-spherical, ADAFs 
are differentially rotating, turbulent, etc.. We suspect that ADAFs lie somewhere between the 
two extremes of thin disks and spherical accretion, but numerical simulations will be necessary 
to determine precisely where. 
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Fig. 1. — Radial profiles of several quantities in our global dynamical models of ADAFs. 
Each panel shows plots for rj from to 1 in steps of 0.2. Other parameters take the values: 
7t = 1.5, a = 0.3,/ = 1, and (3 = 1 at r = r out = 10 5 r 9 . (a) radial velocity (U r ) and total 
sound speed (vt), (b) density, (c) rotation rate in units of the Keplerian rate, (d) ratio of gas to 
turbulent pressure (/3). Note that U r , vt, and p are essentially independent of rj, while 
and, in particular, (3, depend non-trivially on rj. 
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Fig. 2. — (a) (3(r) for the same parameters as in Figure 1, except (3 is set to 15 at the outer 
boundary, (b) f3(r) for rj = 1/2, 74 = 1.5, and several different values of at the outer 
boundary. The solutions roughly converge to a common (3 ~ 5 in the interior (although the 
radial variation of (3 is quite different). 



